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Abstract 



CO \ Parameters of Gasser Leutwyler chiral Lagrangian are proved saturated by 

dynamical quark self energy S(fc 2 ) in a gauge invariant, nonlocal, dynamical 
' quark model. 

o ' 
o 

\ Much of low-energy QCD can be encoded into a series parameters appearing in a chiral 

Lagrangian, expanded to some finite order of low energy expansion. Attempts have been 
made to understand these parameters: it is shown that low lying vector mesons will saturate 
Ph! the parameters ffl. To go beyond phenomenological level, the anomaly contribution was 
taken as the main source of the parameters 0, we call this type investigation the anomaly 
approach, it leads result 

^ ■ N 

8Li = 4L 2 = -2L 3 = 24L 7 = -8L 8 = L 9 = -2L W = (1) 

which are close to experiment result except L-j and Lg which have wrong signs. The defi- 
ciency of this calculation lies in its independence of interaction: if we switch off the strong 
interaction and discuss a system of free quark field with external sources, the anomaly calcu- 
lation can still be performed without any change. Then it seems that ([!]) is not due to strong 
interaction among quarks and gluons, but rather an artificial result. Another type research, 
we call it dynamical approach, mainly consider the dynamical effect in which the main 
source of the parameters is from dynamical quark self energy S(fc 2 ). This approach has 
advantage of maintaining chiral symmetry and momentum dependence of dynamical quark 
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mass, in the mean time avoiding introduce in the theory the hard constituent quark mass 
to cause wrong bad ultraviolet behavior of the theory. But it does not explain why it can 
offer the better numerical result (without wrong sign problem for L 7 and L 8 ) than anomaly 
approach. In fact, anomaly contribution and dynamical quark self energy contribution are 
two independent sources, if the anomaly contribution play role, according to ([I]), it will be 
dominant at all parameters and then there is no room left for S(fc 2 ) to play role to match 
the experiment data, except for L 7 and L 8 . If the anomaly contribution donot play role, it 
must be cancelled in some sense and after the cancellation, we need to show the remaining 
dynamical effect (which may or may not be dominant by dynamical quark self energy) can 
still recover or improve the result ([[]). It is purpose of this work to judge these two possi- 
bilities. We will show the second choice is correct, the cancellation do happen in dynamical 
approach and remanent contribution from dynamical quark self energy S(A; 2 ) can provide 
values for parameters of chiral Lagrangian consistent with experiment data. 

Consider QCD in presence of external scalar, pseudoscalar, vector and axial vector 
sources, 



where ip, are light, heavy and gluon fields respectively. £(ip, i>, ^, A^ is Lagrangian 

of QCD. The chiral Lagrangian relate this generating functional by 



U is pseudo goldstone boson (PGB) field, Sq^U, J] is Gasser and Leutwyler (GL) chiral 
Lagrangian M, 



J(x) = j) (x) + ^ (x)j5 ~ s(x) + ip(x)j 5 , 



The generating functional in Minkovski space is 



Z[J]= / VipV^WV^VA^ e */**[AW>*,*A0-HW] 



(2) 




(3) 




(4) 



+L 2 tr[V M C/ t V„[/]tr[V M C/ t V I '[/] + ^[(V^V^f/) 2 ] 
+L A tr[V fl UW ll U}tr[x ] U + U^ X ] + L 5 tr[V^V ^U( X ] U + U ] X )} 
+L 6 [tr( X t C/ + t/ f x)] 2 + Hto&U - lfx)f + LfctfU^U + X^X^) 
-iLvtT[F*S7»UV v U ] + Fj^VrfVU] + L lo tr[U^F^UF L ^\ 
+# 1 tr[F* F^" + F^ V F L ^} + H 2 ti[x ] x] } + 0{p & ) terms ( 



(5) 
(6) 



Sanomalyt^, J] = S WZW [U, J] + 0(p & ) terms , 
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where 

X (x) = 2B [s(x) + ip(x)} (7) 

and Sy/zw[U, J] is Wess-Zumino-Witten action given in Ref. ||. Up to order 
of p 4 , S^nomaiy J] is completely known, but S noTiaa x[U, J) left fourteen parameters 
F , B , Li, . . . , Lio, Hi, H 2 need to be calculated. To reveal the source of these parameters, 
we improve the conventional dynamical approach by building up a gauge invariant, nonlocal, 
dynamical (GND) quark model. The action in GND model is assumed to be S e g [ip, ip, U, J], 
it relate to our generating functional by 

Z[J]= J VUVipVtf e iSM"P&u,j] _ ( g ) 

The r.h.s. of above equation can be seen as a result of integrating out heavy quark and 
gluon fields and integrate in the PGB field U in (Q). If we further integrate out light quark 
field in above generating functional, we obtain GL result @. So S e s[ip, ip, U, J] can be seen 
as an intermediate stage action to relate fundamental QCD with phenomenological chiral 
Lagrangian. 

Compare (|3|) and (|||), we find GL chiral Lagrangian relate to GND model by 

e is GL [u,j] = J v ^ v ^ e is cS [^,u,j] ^ ^ 

S e g[ip,ip, U, J] is required to be invariant under following local Ul(3) <8> U R (3) chiral trans- 
formations: 

V>(s) -> 1>'(x) = [V R (x)P R + V L (x)P L ]i>(x) 

J(x) - J\x) = [V R (x)P L + V L {x)P R ][J{x) + ip]\Vi{x)P R + Vl(x)P L ] 

U(x) -> U'{x) = V R {x)U{x)Vl{x) . (10) 

Notice that U field has standard decomposition U(x) = Q(x)Q(x) and Q(x) field, under 
transformation fllPp, transform as fl(x) — > Q'(x) = h) (x)Vt(x)Vl(x) = V R {x)Vt(x)h(x) with 
h(x) depend on V R , Vl and Q, represent an induced hidden local U (3) symmetry to keep 
transformed Q be a representative element at coset class. 

To implement local chiral symmetry explicitly, we take a special local chiral transforma- 
tion V R (x) = Q^x), Vl(x) = Q(x), the corresponding hidden symmetry transformation is 
h{x) = 1, 

<Mx) = [Q\x)P R + Q(x)P L ]tfj(x) (11) 
J n (x) = p(x)P R + n\x)P L ] [J{x) + ip] p(x)P R + &(x)P L ] 

= sn(x) + ipn(x)l5 + f n(x) + </i n{x)^ (12) 
U a (x) = 1 . 
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On rotated basis, we can rewrite (|9T) as 



JVtpVtpe i -f d4x ' p{xWx+Jix)]i;ix) J 



(13) 



where N' = fV^V^e i S^^' +J ^^ = Det[ip x + J(x)}. In the last equality, we have 
taken chiral rotation (|TT|) for functional integration measure both in numerator and denom- 
inator. The possible anomalies caused by this rotation are cancelled between numerator and 
denominator. Since we are only interested in U dependence of the theory, pure source terms 
N' is irrelevant and therefore can be treated as a normalization factor. 
Result ( |13"D tells us that S c s should has following structure 



S eS [ipn,ip n ,l,JQ) = J d*xifj n (x)[ip x + J n (x)]ip n (x) + S int [ip n , ^n, h M ( 14 ) 

where Si n t is interaction part caused by color gauge interaction. If we switch off color 
gauge interaction which means we are dealing with free fermion fields, there will be no 
effective Lagrangian (Sql = 0). Si n t should include those fermion self interaction terms 
caused by integrate out gluon and heavy quark fields in underlying QCD and integrate in 
local goldstone boson fields U. Among these, the most important effect related to chiral 
symmetry at low energy region is spontaneous chiral symmetry breaking (SCSB) which 
require quark has a nontrivial momentum dependent self energy S(fc 2 ), its effects can be 
introduced into the theory by adding in S- int a self energy term 

d 4 x^ n (xMdl)M^) ■ (15) 

Just this term itself is not enough, since it is not invariant under local chiral symmetry trans- 
formations. To make it invariant, in conventional dynamical approach 0, a non-integratable 
face factor is introduced into theory which cause very complex formulae and authors in 
even donot put their analytical result in their papers. We donot use non-integratable face 
factor, instead we note that local chiral symmetry transformation on rotated variable is 

Jn(x) -> J' n (x) = rf(x)[J n (x)+ip x ]h(x) ■ (16) 

Original local chiral symmetry now is realized as a hidden local symmetry. Once the theory 
is constructed to be invariant under this hidden symmetry, it is invariant under original 
local chiral symmetry. Since interaction part in S e ff should be invariant on local chiral 
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symmetry, we need at least to generalize self energy term (|15|) to be invariant on hidden local 
symmetry (Jl6|). To achieve this, we change the ordinary derivative d£ to hidden symmetry 
covariant derivative = — Wq(x) (the overline on is to denote the difference with 
covariant derivative appeared in (|^)). ([IB]) tells us Vq(x) transform as Vq(x) — > Vq (x) = 
hi (x)vq(x)Ii(x) + ih)(x)[d^h{x)} which lead — > V^' = W{x)V^,h{x). The modified chiral 
invariant interaction action now is 

s int [^ n , i, M = -J d 4 % M*Mv 2 x )Mx) • (17) 

This action is not the complete part of interaction, but it is the minimal part of interaction 
which respect local chiral symmetry with dynamical quark and SCSB. If we take the idea 
of dynamical perturbation originally from Pagel-Stokar || and developed in Ref. 0, in 
which at the leading order of the expansion, all perturbative effects are ignored and only 
nonperturbative effect considered in the theory is that from quark self energy E(/c 2 ). ( |l7|) in 
this sense can be seen as a result of leading order expansion from dynamical perturbation. 

In GND model, quark fields dependence is bilinear and can be exactly integrated out, 
the result GL Lagrangian from ( p3| ) is 

S GL [U, J] « S G nd[U, J] = -tTr ln[ip + J Q - S(V 2 )] + zTr ln[ip + J a ] . (18) 

Use the Schwinger proper time formulation developed in ||, we can compute the E(V ) 
dependent determinant in (|i~8"D. The result is 



-ili ln[ip + Jn- S(V 2 )] (19) 
= J d 4 xtr f B FqSq + Cia 2 ^ + C 2 [d^a^} 2 + C z (d^a v n - d u a^)(d^a n ,u - d v aa tli ) 



+Cupnd^a n 



+ 0(p % ) + imaginary terms 



where tr/ is trace for flavor indices. Covariant derivative for function / and Vq are defined 

as 

#/ = &>f - ivjtf + zK V/T = d^ n - - z« + « , ■ (20) 

The £ dependence for coefficients appeared in ( |19"D are 

F*B = 4 1 dkZ k X k (21) 

(22) 



d = 2 / dk 
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Xk , k w2 X fc 



C 2 = -2fdk[- 2A k X\ + 2^ - + '-Y' 2 ^ - ^£' fc 2 X, 

C 3 = - J dk 



X 2 

- 2B k X k + - ^ a 4 + Y k A 2 " ~ 2 kk 



dk 



,4£ 4 2k 2 Y\ k\ 4 6X| 
'~3 3~ + 18 )(6Xfc "^ 



^ y t 2 
+ ^-§) + (-4^ + ^)(-2X| 



2X? Xt. . Xi. 

+ -^-71—^ + ^? 



C, = / (ifc 



A 2 A 4 ' A 2 

-4£ 4 2A; 2 S 2 k A 

( H H 

V 3 3 18 

X k X k 2 

"a^+a 2 " Afc 



/:■ -V- , 2: )(6X^ i ^ feN iv2/ ov:! "" x 



A 2 



A 4 



9 V 2 



C 6 = 2Jdk 
C 7 = 2 [ dk 



(3£ 2 + 2A;%£' fc )X 2 + [-2E 2 - A; 2 (l + 2£ fe £' fc )]^f 

Xi. 



^ A 2 

X 2 



(S 2 +2A; 2 S^)X fe 2 -A; 2 (l + 2S fc S / A 
C 8 = 4/ ^[(-4E 3 fe + A; 2 £ fc )Xf + (4£ 3 - k 2 Y k )^ - (2S 3 fc - \k 2 Y k )^ + 3S fc ^f - 3£ fc X 2 



Xu 



dk 



,l -k 2 Y' k Y'l + ±E k TT k )X k + (C fe - £> fc )§ - (C k - D k )X 2 - 2E k Xl 



+2E M-E k ^ 
A 2 A 4 



iCio = 4: J dk 
C11 = -4jdk 
where 

y ^ = ^ c y 



2F fe Xf + 2F k § - F fc ^ + ^Y' 2 ^ - ^£' 2 X 2 
-(^ + ^)|| + (S fc + iA; 2 ^)X 2 



(23) 



d 4 fc fc 2 -s 2 (-fc 2 ) 

rfA; = iX r / — e ^ 



(2tt)< 



B k 

C k 
D k 



n-k 2 ) 
1 



X*. 



A; 2 - S 2 (-/c 2 ) 



^ 2 S fc S^l + 2S fc S' fc ) + is 2 (l + 2Y k Y' k ) - ^ 2 S 2 (S^ 2 + + h\Y' 2 + £ fe E£) 

^ 2 S fc S^(l + 2Y k Y' k ) + is 2 (l + 2Y k Y' k ) - ^ 2 S 2 (S^ 2 + £ fc ££) + + E*E£) 

+^ 2 (l + 2S fc S' fc ) 

1 1 1 

- - -Y,Y' - -k 2 Y' 2 

3 3 LfeLfe 2 k Lfe 

1 -,2^/2 , 1 l2 



2 * 2 E* + ^£^(1 + 2Y k Y' k ) + ^'^'(l + 2Y k Y' k ) + ^(S^ 2 + E fc E£) 
+ 1 -k 2 Y k Y' k (Y' 2 + Y k Yl) 
E k = -h 2 Y k Y' k {\ + 2S fc S^) 2 - h 4 Y' 2 (l + 2£ fc ^) 2 
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Fk = ~^kK + h\^' k r - hi + 2 -Y? h Y,' k - h^lVZ + E fc E2) + ^ 4 (£' fc 2 + E fc E£) 
+ifc 2 (l + 2£ fc £' fc ) - h 2 . 
The result for C\ = _F 2 in ( p2|) is just the well known Pagel-Stokar formula 0, if we take 
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momentum cutoff A be infinity. The part of S(V ) independent quark determinant in fll9|) 

is just the result of anomaly approach with a total minus sign, we can get result of anomaly 

approach by taking limit of = const — > (Note due to possible infrared divergence, limit 

of E — > must be taken after the momentum integration). The nonzero coefficients Ci for 

the case of infinite momentum cutoff A is 

„ N c „ N c „ £ 2 v e-o AT C „ S 2 

ft — 24^. ft — 48^(^+7 + D, ft — ^(^+7 + 4), 

~ #c S 2 £->o iV E^O iV 

„ s^o N c , S 2 „ s^o N c /f S 2 

ft — - 4g^(i»X5 + 7) ' * Cl ° — 12^ ( ^ + 7 + 2) ' 



The pure imaginary terms in (|19|) are completely known at phenomenological level and 
its calculation in terms of E is already performed in Ref. |J and proved exactly recover the 
Witten's result ||, we donot explicitly write down their detail structures. 

With help of ([12]), ([18]) and (0) lead relation, 

1 C5 Co .C10 r C4 — 2C5 + 6C9 — 3iCio 
Li = — Lo = h 2 , L s = , 

2 32 _16 32 _ _16 

L 4 = 0, L 5 = —7— , Lq = 0, L 7 



16£ 48 485 

r _ft j ft _ £7 C11 . -4C 9 + zgip 
8 ~ 16 + 165 2 16£? 2 + 16B ' 9 ~ 8 

_ -C 3 + C 9 _ C 3 + C 9 _ C 2 C% C7 fti m ^ 

Lw ~ 2 ' Fl " ~ i - ' ^"T + sif + sif - ^- (24) 

where Cj = Cj — lims-*oCi % = 1,2,.. .,11. Term — lims-+oCj is of special interest, since 
it relate to anomaly result mentioned in the beginning of this paper. In fact, in anomaly 
approach, the effective action is iTr ln[if) + Jq] 0, which is just the result of fl5|) with 
(f2~4|) by taking C\ values at E = and revert all signs. One can easily check this reproduce 
result (0) in which all ultraviolet divergence are cancelled each other for Lj parameters. The 
interpretation of this result is that from (|i~8D , the contribution of anomaly play no role in the 
final result , it is completely cancelled by dynamical quark self energy dependent part, only 
the remainder after cancellation play role in the final parameters Li . 

In (^), parameter B needs special treatment. Since with help of (PH), we find FqB is 
generally divergent. To renormalize this condensate, we note that F and m^ipip)^ (m\ is 
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bare current quark mass) is renormalization invariant or more general, the x field denned 
in (0) is renormalization invariant, i.e. x( x ) = 2B [s(x) + ip(x)] = 2B r [s r (x) + ip r (x)] 
with B r ,s r ,p r are renormalized B Q ,s,p. Correspondingly, renormalized quark condensate 
(ipip)r is (ipip)r = —NfF^Br. So replacing the scalar and peudoscalar sources and B with 
renormalized ones donot change value of x held. With renormalized sources, we can replace 
B in ( p4|) with B r . In this paper, we donot directly calculate and use B , instead we 
calculate and use B r . The renormalization point is chosen to be at scale of lGeV. 

Now, once the quark self energy S(fc 2 ) was input into the formulae, we can get all 
parameters in Gasser-Leutwyler chiral Lagrangian. In conventional dynamical approach || , 
the ignorance of E(/c 2 ) is parametrized by following ansatz 

which satisfy E(m 2 ) = m and shares qualitative similarities with solutions of improved ladder 
SD equation. It is finite, positive, monotonically decreasing functions with 1/k 2 behavior 
at large q 2 and S'(0) < 0. The constituent quark mass m is determined for each choice of 
A for F = 93MeV from (f^). For A = 1,2,3,4, we obtain m = 379, 350, 331, 3l7MeV 
respectively. With ansatz (^), the result parameters are listed in TABLE.|. We see that 
the wrong sign problem for L7 and L§ in anomaly calculation is corrected now and result 
parameters are roughly consistent with experiment data. 

To further trace the relation of GND model with underlying theory QCD. Note that 
the parameters in the chiral Lagrangian are recently expressed in terms of QCD Green's 
functions |I| and for quark two point Green's function $q P (x,?/), at large N c limit, flO 



gives equation, 

i) n+1 (N c g 2 ) n 



[1$ + i^f x + J n + Z] ap {x, y) + ]T / d i x l ■ ■ ■ d*x n d*x[ ■ ■ ■ d 4 x' r 

n=l J 

fa y,xx,x[,---, x n , x' n )<S>£ pl ( Xl ,aQ ■ ■ • (x n , x' n ) = 0. (26) 

Now consider the coincidence limit of two point quark Green's function in QCD in 
presence of external sources, 

^<0|TCW^fa)|0>QCD-~* 6lnZ[J] 



_ 1 JVU S j?%!%J x) e iSGM _ JVU $H)( fe €)[[/,J](x,x)e^^ J ] 
" N c JVU e iS G^ u > J ] ~ JVU e i5 G L [tV] 

where we have used (H) and c^MOO [U, J](x, x) is 



(27) 



S 



[ ' R ' j " N c 8J(^)(bO( x ) N c 5J^m)( x )- 



Use (p~8|) , the rotated $ become 

^ vm) [U, J]{x,x) = Tr 

which imply 



5[Jo - S(V 2 ) 



Jn-S(V 2 ) 54 ar?)W) (x) 



(29) 



$ n (x, y) = -ij d*z \ip + J n - J:(V 2 )]-\x, z) [5(z -y) + A(z, y)} (30) 



where A(z,y) relate to <5£(V )/SJ. Compare to (p6|), we find, present choice of Si nt is 
equivalent to take following approximation 



l + Ar x [E(V 2 )-A(^ + J Q )] 



ap 



{x,y) 



n=l ' 



-"-/',' , „ i .... V /,/>, ...,/' r ,/ »,-'... J4 X ' V *)" + ( N c9 

n! 



(x, y) + / <i 4 xi ■ ■ • c? 4 x n c? 4 a;i • ■ ■ <i 4 



xc;^. ■ • • . ^. < (xi, «i) ■ ■ ■ o . (si) 



PP\-"Pn 

yCr<7 1 ■■■CT, 1 



If we further drop correlation functions G ° (x, ?/, Xi, ■ • ■ , x n , x' n ) with n > 1 and ignore 
the external sources in above equation (A therefore must be ignored), as mentioned in Ref. 
IIJ, pi| ) then is just Schwinger- Dyson equation for quark propagator. We only consider the 
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term of quark self energy with argument of V which is the minimal generalization from pure 
X(<9 2 ) term to source dependent terms satisfying local chiral symmetry. Include in source 
terms, just self energy term in l.h.s. of ( |3~TD is not enough to match all contributions of r.h.s. 
of equation, but if underlying QCD can provide correct predictions for parameters in chiral 
Lagrangian, the fact of S(fc 2 ) dominance in the parameters of chiral Lagrangian imply that 
the terms we dropped at this work should not play so important role, we will leave these 
terms in addition to self energy term to balance equation ([H]) in future investigations. 

In conclusion, all 12 parameters in p 4 order SU(3) chiral Lagrangian are explicitly ex- 
pressed in terms of functions of quark self energy S(/c 2 ). We have shown that the original 
result given from anomaly approach is completely cancelled. Take suitable quark self energy 
ansatz to perform numerical calculation, we find after cancellation of anomaly contribution, 
the dynamical quark self energy do can provide parameter values roughly consistent with 
experimental data. 
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TABLES 

TABLE I. Values, multiplied by 10 3 , for the parameters of the order p chiral Lagrangian 



calculated in GND model with quark self energy determined by ansatz (25). anomaly: anomaly 
approach result; expt: experimental values. 





Li 


L 2 


L 3 


u 


L 5 


L e 


L 7 


^8 


L 9 


Lio 


A=l 


0.927 


1.85 


-6.55 





1.79 





-0.570 


1.56 


3.79 


-5.18 


A=2 


0.771 


1.54 


-5.53 





1.63 





-0.501 


1.39 


2.60 


-3.91 


A=3 


0.708 


1.42 


-5.09 





1.51 





-0.449 


1.26 


2.13 


-3.38 


A=4 


0.674 


1.35 


-4.85 





1.41 





-0.413 


1.17 


1.90 


-3.08 


anomaly 


0.792 


1.58 


-3.17 











0.263 


-0.792 


6.33 


-3.17 


expt 


0.9±0.3 


1.7±0.7 


-4.4±2.5 


0±0.5 


2.2±0.5 


0±0.3 


-0.4±0.15 


1.1±0.3 


7.4±0.7 


-6.0L0.7 
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